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ABSTRACT 

The stability of metal-free very massive stars (Z = 0; M = 120-500Mo) is analyzed and compared with metal- 
enriched stars. Such zero-metal stars are unstable to nuclear-powered radial pulsations on the main sequence, 
but the growth time scale for these instabilities is much longer than for their metal-rich counterparts. Since they 
stabilize quickly after evolving off the ZAMS, the pulsation may not have sufficient time to drive appreciable mass 
loss in Z = stars. For reasonable assumptions regarding the efficiency of converting pulsational energy into mass 
loss, we find that, even for the larger masses considered, the star may die without losing a large fraction of its 
mass. We find a transition between the e- and K-mechanisms for pulsational instability at Z ^ 2 x 10"'* - 2 x lO"-'. 
For the most metal-rich stars, the ^-mechanism yields much shorter e-folding times, indicating the presence of a 
strong instability. We thus stress the fundamental difference of the stability and late stages of evolution between 
very massive stars born in the early universe and those that might be born today. 

Subject headings: stars: very massive — stars: pulsations — stars: mass loss — stars: zero metallicity 



1. INTRODUCTION 

The formation and nature of the first generation of so-called 
Population III (Pop III) stars have been studied and speculated 
about for over thirty years (Schwarzschild & Spitzer 1953; 
Yoneyama 1972; Hartquist & Cameron 1977; Palla, Salpeter 
& Stabler 1983; Silk 1983; and many others), but the extent 
to which they differed from present day stars in ways other 
than composition is still debated. Recently, three-dimensional 
cosmological simulations have reached sufficient resolution on 
small scales to begin to address the star formation problem (Os- 
triker & Gnedin 1996; Abel et al. 1998). While considerable 
uncertainty remains regarding the continued evolution of the 
dense knots they find in their calculations, these simulations do 
not exclude the formation of a first generation of quite mas- 
sive stars - 100- 1000 Mq (e.g., Larson 1999; Abel, Bryan, & 
Norman 2000). Using different numerical methods to analyze 
the fragmentation of primordial clouds, Bromm et al. (1999) 
reached similar conclusions. Without further study, especially 
including the complexity of radiation transport and molecu- 
lar opacities, it may be premature to conclude that all Pop III 
stars were so massive, but there certainly exists adequate mo- 
tivation to examine the properties of such stars. In particular, 
if such massive Pop III stars were ever born, would they re- 
tain their large masses until death? That is, would they die as 
pair-instability supernovae, with possible unique signatures of 
nucleosynthesis and black hole formation, or would they lose 
most of their mass and die much the same as present day stars? 

It is well known (Ledoux 1941; Schwarzschild & Harm 
1959) that, above a critical mass, main sequence stars are vi- 
brationally unstable due to the destabilizing effect of nuclear re- 
actions in their central regions. For stars of Pop I composition, 
the critical mass for this instability is ^ lOOM© (Schwarzschild 
& Harm 1959). According to non-linear calculations, such an 
instability leads to mass loss rather than catastrophic disrup- 
tion (Appenzeller 1970a; Talbot 1971ab; Papaloizou 1973ab). 



Consequently, it has been thought that stars more massive than 
about 100 Mq would not survive for long before losing a sub- 
stantial fraction of their mass. Although the stability of very 
massive stars has well been explored since the work of Ledoux 
(1941), somewhat surprisingly, all these analyses have been 
devoted to solar-like compositions or slightly metal-poor stars 
(Schwarzschild & Hiii-m 1959; Stothers & Simon 1968; Aizen- 
man et al. 1975; Stothers 1992; Glatzel & Kiriakidis 1993; Kiri- 
akidis et al. 1993), or to helium-core objects such as Wolf-Rayet 
stars (Maeder 1985; Cox & Cahn 1988). Since the temperature 
dependence of hydrogen burning and the surface gravitational 
potential of low metal stars are quite different, it is important to 
learn whether their stability criteria may also vary. 

We have thus performed, to our knowledge, the first stabil- 
ity analysis of zero metallicity stars with 100 < M/M© < 500 
and have also studied the stability of stars in this mass range 
for variable initial metal mass fraction < Z < 0.02. The stel- 
lar models are presented in Sect. 2. The results of the linear 
stability analysis are given in Sect. 3, mass loss estimates are 
presented in Sect. 4, and our conclusions follow in Sect. 5. 

2. STELLAR EVOLUTIONARY MODELS 

The evolution of four stars of 120, 200, 300 and 500 M© was 
followed from the Zero Age Main Sequence (ZAMS) to the end 
of core helium burning using an implicit Lagrangian hydrody- 
namics code adapted to stellar modeling (Langer et al. 1988, 
Heger 1998; Heger et al. 2000a). Details of the code can be 
found elsewhere (Heger & Langer 2000), but to summarize the 
main physics, the opacity is taken from Iglesias and Rogers 
(1996) and the principal nuclear reaction rates from Caughlan 
& Fowler (1988) with the rate for '^C(a,7)"'0 multiplied by a 
constant 1.7 (see also Buchmann 1996). Mass loss is not taken 
into account, except for what is driven by the pulsations studied 
here. For zero metallicity stars, radiative mass loss is probably 
negligible (Kudritzki 2000). 
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Fig. 1 . — Theoretical H-R diagram for stars with Z = and different masses 
assuming no mass loss. The time-scales of evolution are indicated by marks 
on the curves. 



the hydrogen-burning shell becomes more active with increas- 
ing stellar mass. 

In order to analyze the pulsational stability as a function of 
metallicity, we repeated our calculations for a wide range of 
Z from Z = 2 X 10"'' to solar, Z = 0.02, for the 120 Mq, and 
for selected values of Z > forM > 120Mq. The results are 
displayed in Fig. 2 for a 120 M0 star. Since most of the stabil- 
ity analyses mentioned in Sect. 1 (Ledoux 1941; Schwarzschild 
& Harm 1959; Appenzeller 1970a; and others) are based upon 
models with solar metallicity, but computed using only electron 
scattering opacity, we also calculated, for comparison, models 
with such an opacity (Fig. 2). 
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Fig. 2. — Theoretical HR diagram for 120Mq stars with different Z from 
hydrogen ignition till I % of hydrogen left in the center of the star - assuming 
no mass loss. [Z] corresponds to log(Z/Z0) 



As was first noted by Ezer and Cameron (1971), massive 
stars with no metals do not obtain their energy on the main se- 
quence from the pp-cycle as one might naively expect. Instead 
the star contracts to >10**K, burns a trace of helium to carbon, 
and then expands slightly to burn hydrogen by the CNO cy- 
cle. A CNO mass fraction of just a few times 10"^ is sufficient 
to catalyze even the enormous energy generation of such stars. 
However, owing to this small CNO mass fraction, the star re- 
tains a central temperature of ~ 10*^ K throughout central hydro- 
gen burning. As Ezer & Cameron (1971) also noted, these stars 
are thus much more compact than their metal-enriched counter- 
parts and their radiation, while on the main sequence, is in the 
ultraviolet (Table 1). 

Figure 1 shows the theoretical tracks of our Z = stars in 
a Hertzsprung-Russell diagram (HRD). All begin to contract 
dynamically after core helium burning and the star ceases to 
evolve in the HRD due to shorter evolutionary time-scale of the 
core compared to the thermal time-scale of the envelope. Note 
however that, forM>300MQ, even metal-free stars evolve to- 
ward Teff < lO'^K and eventually become red supergiants, as 



3. LINEAR STABILITY ANALYSIS 

A linear non-adiabatic stability analysis of all models was 
carried out using a radial pulsation code originally developed 
by Umin Lee (1985; see also Heger et al. 1997; Baraffe et 
al. 1998; Alibert et al. 1999 for recent applications). The main 
uncertainty in these calculations is the treatment of convection, 
which is treated in the "frozen-in" approximation, i.e., the per- 
turbation, SFcom, of the convective flux is neglected in the lin- 
earized energy equation. Although this approximation is an un- 
avoidable consequence of our poor knowledge of the interaction 
between convection and pulsation, it should not be crucial. As 
we shall see, the relevant regions for pulsation, which are lo- 
cated in the deep interior of the stars, are fully convective and 
can be considered as adiabatic. In the outer regions, which may 
also be important for the pulsation, the convective flux is not 
the dominant energy transport. 

The stability of the fundamental mode and the overtones (up 
to the fifth overtone) was analyzed starting from ZAMS models 
up to the end of core helium burning for Z = stars. 
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Fig. 3. — (a) Differential work dW/dM, in ai'bitrary units, as a function 
of tlie temperature (T in K) in the interior structure of a 300 Mq with Z = 0, 
on the ZAMS. The Rosseland mean opacity k (in cm" g" ) is indicated by the 
dash-dotted line. The dashed line coiTesponds to the nuclear energy generation 
e (in units of 5 X lO' erg g" 
mass. 
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Fig. 4. — Stability coefficient, K, as a function of time. Positive values in- 
dicate growth of the pulsation amplitude while negative values (gray region) 
correspond to damped oscillations. Shown are 120 Mq (solid line), 200 Mq 
(dashed), 300 Mq (dash-dotted), and 500 Mq (dash-triple-dotted) stars of 
Z = (assuming no mass loss). The filled circles indicate the time where these 
stars become pulsationally stable. 



3.1. Metal-free stars 
3.1.1. ZAMS models 



Table 1 gives, for metal-free ZAMS stars, the effective tem- 
perature, Teff, stellar radius, R, central temperature, T^, pulsation 
period, P, and the pulsational e-folding time, t^ (Cox 1980). 
Adopting a time dependence for the eigenfunctions of the form 
exp(icrf ) exp(Kt), where a = Itt/P is the eigen-frequency and 
K, the stability coefficient, the e-folding time for the growth 
of the instability is given by ra = 1 /K. For all masses studied, 
the ZAMS models are unstable and oscillate in the fundamental 
mode, whereas the overtones are stable. The excitation mech- 
anism is related to the temperature and density dependence of 
the nuclear energy generation rate near the center of the star, the 
so-called e-mechanism. As was shown by Ledoux (1941), this 
driving mechanism is efficient because of the large radiation 
pressure in such massive stars. The driving and damping zones 
are indicated by correspondingly positive and negative values 
of dW/dM, where W is the work integral (Cox 1980; Unno et 
al. 1989). This is illustrated in Fig. 3 for the 300 M© star. Note 
that the outer damping regions (see Fig. 3a) give a negligible 
contribution to the work integral due to the very small amount 
of mass involved (see Fig. 3b). The same behavior of the work 
integral is also found for the other masses studied. 



3.1.2. Effects of evolution 

It is known from previous studies of solar metallicity stars 
that high mass stars eventually stabilize as they evolve away 
from the ZAMS (Schwarzschild & Harm 1959; Stothers and 
Simon 1968; Aizenman et al. 1975; Maeder 1985). Stabiliza- 
tion comes about because of the increase of the mean molecular 
weight in the inner regions which increases the central conden- 
sation (see, e.g., Maeder 1985 for details). The same properties 
are also found here for metal-free stars. As they evolve away 
from the ZAMS, their e-folding time increases, indicating a de- 
crease of the driving term due to nuclear energy compared to 
the damping term from heat leakage in the envelope. This is a 
direct consequence of the increase of the central density relative 
to the mean density. The stars become stable (negative K) on 
the main sequence after a decrease of the central hydrogen mass 
fraction. This is illustrated in Fig. 4, which displays the stabil- 
ity coefficient, K, as a function of time for different masses, and 
in Fig. 5, which shows the time integration of /T during core hy- 
drogen burning. In Fig. 4, the larger value of K at early times 
and its rapid variation with time are due to the adjustment of the 
stellar structure when the initial CNO is produced (see §2). The 
results are also summarized for the different models in Table 2 
which gives the number of e-foldings, A^ax, by which the am- 
plitude of an initial perturbation would grow in the framework 
of linear stability analysis, as derived from integration of the 
stability coefficient with time (see also Fig. 5). 
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Fig. 5. — Time integration of tlie stability coefficient, K. Sliown is tlie in- 
tegral ^'^K(t')ik' for 120 Mq (solid line), 200 Mq (dashed), 300 Mq (dash- 
dotted), and 500 Mq (dash-triple-dotted) stars of Z = (assuming no mass 
loss). This gives the number of e-folds by which an initial pulsation would 
grow till that time. The filled circles indicate the time where the stars become 
pulsationally stable. After that the pulsation amplitude is damped quickly. The 
slope of the curves is shown in Fig. 4. 



After reaching pulsational stability during central hydrogen 
burning, all stars studied remained stable for the rest of cen- 
tral hydrogen burning and most of central helium burning as 
well. Stars of < 200 M© remained stable throughout core he- 
lium burning, but stars of initial mass >300Mq became unsta- 
ble again towards the end of helium burning. This resurgence of 
instability happens for a central helium mass fraction Fc = 0. 1 in 
the 300 M0 star and Y^ = 0.2 for the 500 M© star, correspond- 
ing to 30 and 50 kyr before the end of central helium burning, 
respectively. The unstable mode is the fundamental mode, char- 
acterized by P = 14. 2d andTd = 0.13yrforthe 300 M0 star and 
P = 40d and t^ = 0.16yr for 500 M© star. The pulsations in 
these blue supergiants, with T^ff ~ 10"* K, exist only in a surface 
layer that contains very little mass. The 500 Mq star, however, 
expands further and becomes a red supergiant with an extended 
convective envelope that is characterized by increasing period 
and decreasing e-folding time. 

The instabilities during core helium burning are character- 
ized by small values of rj compared to the pulsation period. 
This indicates the presence of a strong instability that can de- 
velop rapidly and reach large amplitudes. Inspection of the 
work integral shows that the driving zone is located in the hy- 
drogen/helium ionization zone and the excitation mechanism is 
the standard K-mechanism. 

Note also that the helium cores would quickly become unsta- 
ble to the e-instability, just like their higher metallicity counter- 
parts, if the hydrogen envelope were removed, e.g., to a binary 
companion. 



3.2. Effect of Z 

Table 1 gives the pulsational properties of ZAMS models as 
a function of Z. All models analyzed for Z > are vibra- 



tionally unstable. The fundamental mode is the pulsation mode 
and the overtones are stable. Only for 120 M© with Z > 0.01, 
are both the fundamental (F) and the first overtone (IH) unsta- 
ble, with an e-folding time for IH much shorter than for the F 
mode. Table 1 also gives the results for solar metallicity mod- 
els calculated with pure electron scattering opacity. For such 
models, our results are in good agreement with previous results 
based on pure electron scattering or on simple opacity laws. 
Schwarzschild & Harm (1959) find for 121 Mq, P = 0.38d and 
Td = 1800 yr, and for 218 Mq, P = 0.54d and t^ = 930yr. Aizen- 
man et al. (1975) obtain for 1 IOMq, P = 0.37d and Td = 644 yr. 

Inspection of Table 1 shows that for a given mass, the period 
decreases with metallicity. This is a consequence of the smaller 
radius and more compact structure of metal-poor stars when 
they reach thermal equilibrium on the ZAMS, a well known ef- 
fect of Z on evolutionary models (e.g., Ezer & Cameron 1971; 
Baraffe & El Eid 1991; Meynet et al 1994). The per iod, P , 
is related to the dynamical time scale of the star, P oc ^JWjM. 
Another important result of Table 1 is the increase of Td as Z de- 
creases from Z = 0.02 to 0. This indicates that metal-free stars, 
although unstable, require a much longer time for the growth of 
the oscillation amplitudes than solar metallicity stars. 

We also note a significant difference between Td obtained 
for Z = 0.02 models calculated using the most recent opaci- 
ties, and the electron scattering cases. This highlights an im- 
portant contribution of opacity to the instability of very mas- 
sive stars with solar metallicity. Indeed, the driving zones in 
the Z = 0.02 ZAMS models with M > 2OOM0 are located at 
logT ~ 5. 2-5.3 K and are associated with a sharp peak of the 
opacity in this region due to the contribution of heavy elements 
(Iglesias & Rogers 1996). For I2OM0 with Z > 0.01, our 
results are in agreement with the analysis of Glatzel & Kiri- 
akidis (1993): the driving for the first overtone is due to the 
K-mechanism whereas for the fundamental mode, both the e- 
and the k- mechanisms contribute to the driving. 

As Z decreases, we find a transition between nuclear-driven 
pulsation and opacity-driven pulsation (i.e., the e- and k- 
mechanisms). The e-mechanism is dominant for Z<2xl0~^ 
and the K-mechanism dominates for Z>2x 10"^, for the mass 
range of interest. The increase of Td as Z decreases for 
Z<2xl0~^ is partly due to the different temperature depen- 
dence of the nuclear energy generation from the CNO cycle. As 
shown in Table 1 , central hydrogen burning proceeds at higher 
T for decreasing Z. 

Locally, the temperature dependence of the specific nuclear 
energy generation rate, tnuc, depends on temperature as Cnuc oc 
T", where jy := dlncnuc/dlnT. Since the temperature sensitiv- 
ity, ly, of hydrogen burning by the CNO-cycle decreases with in- 
creasing temperature (e.g., Kippenhahn & Weigert, 1992, Fig- 
ure 18.8) the metal -rich stars have a stronger sensitivity than 
the metal poor stars. We find i^ ~ 10 for Z = and i^ increases 
with Z up to -- 14 for Z = 0.02 in the 120 M© star Analytically 
ly can be computed from ly = (t-2)/3 where 



T = 4.248 



ZjZj 



A^Ao 



-/T, 



1/3 



(1) 



-A1+A2 

Al, A2, Zi, Z2 are mass number and charge number of the two 
input nuclei, and Tg is the temperature in 10^ K (Clayton 1968). 
For "'N(p,7), i^ = 14 at Tg = 0.04 and ly = 10 at T<) = 0.1, in 
agreement with the values obtained above. 

Since the driving of pulsations due to the e-mechanism de- 
pends directly on ly with W oc z^, it is larger in metal-rich than 
in metal-free stars, contributing to the decrease of Td = 1 /K oc 
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1/W as Z increases. Another contribution to t^ comes from the 
total pulsation energy fipuis oc 1 /P^, since K oc I //ipuis oc P^ 
(see next section Eqs. (3) and (4)) and P increases with Z, as 
mentioned above. As long as the contribution of metals to the 
opacity is negligible, this explains why t^ = 1/K oc 1/P^W de- 
creases as Z increases. Once the K-mechanism operates in the 
envelope, it yields a dramatic increase of W and significantly 
smaller rj, by up to six orders of magnitude (see Table 1). 

4. MASS LOSS RATES 

In order to derive a realistic mass loss rate, non-linear cal- 
culations would be required to determine the pulsation ampli- 
tudes and dissipation due to shock wave formation. Non-linear 
calculations have been performed for the solar metallicity case 
by various authors (Appenzeller 1970ab; Ziebarth 1970; Pa- 
paloizou 1973a; Talbot 1971ab) who, unfortunately, give dif- 
ferent results. Appenzeller (1970ab) finds an upper limit M = 
5xlO"^M0yr-' for a 130 Mq star andM = 5xlO"4M0yr-' for 
a 270 M0 star, whereas Papaloizou (1973b) finds stronger dis- 
sipation effects and mass loss rates less than 10~*'MQyr"' for 
similar masses. 

Since the dissipation effects in Appenzeller (1970ab) may 
be underestimated compared to other studies, yielding larger 
mass loss rates, we used his results to derive an upper limit 
for M for the Z = models. Appenzeller found that when the 
velocity amplitude at the surface approached the sound speed, 
mass loss sets in. An upper limit of the mass loss rate came 
from assuming that the entire rate of gain of pulsation en- 
ergy, Lpuis, was used to accelerate the matter to escape velocity. 



.JlGM/R. Thus one has: 



M 



-'puis 



(2) 



Lpuis can be estimated from the total pulsation energy fipuis, 
averaged over one period, and the stability coefficient, K (see 
Cox 1980), with 

ipuis = 2/riipuis , (3) 

where 

£puis=2fyj j S'^dM. (4) 

Since the relative pulsation amplitude, 5r/5rsnrS, is given by 
our linear stability analysis, dr^ud and thus iSpuk are estimated 
by assuming, like Appenzeller (1970ab), that the surface veloc- 
ity amplitude reaches the sound speed, Csound, i-c: 



2^. 

p ^^surf ~ *-sound • 



(5) 



For the zero metallicity stars from 120 to 500 M0, Csound is 
typically (1 - 1. 6) x 10^ cms"'. The resulting energy, fipuis, and 
mass loss rates for the Z = models on the ZAMS are given 
in Table 3. Results are also given for the 120 Mq star with 
Z = 0.02 and pure electron scattering opacity. Our upper limit 
for M in this case is in good agreement with Appenzeller's value 
of 5 X 10"^ Mq yr"' for similar masses. A crude upper limit for 
the total mass lost /SM during hydrogen burning is derived as- 
suming that mass loss proceeds at constant rate as long as the 
star remains unstable, i.e., during the time Tstabie (see Table 3). 
In any case, AM does not exceed 5 % of the total mass for zero- 
metal stars. Note that this estimate neglects the time required 
for the pulsation amplitude to grow to the value Sr^arf character- 
izing the beginning of mass loss (Eq. (5)). Given the non negli- 
gible value of Td compared to TstaWe (see Table 1), this time may 



be significant for the 120 Mq star. In fact, this star may never 
encounter pulsationally driven mass loss. Moreover, since the 
stars stabilize as soon as they evolve away from the ZAMS (see 
Sect. 3.1.2) with decreasing K (see Fig. 4), one can expect de- 
creasing M along the main sequence, and thus even smaller val- 
ues for /SM. We thus find that even for 500 Mq star, pulsations 
do not lead to significant mass loss on the main sequence. 

For post-main sequence evolution, as found in Sect. 3.1.2, 
only stars with initial masses of >3OOM0 become unstable 
again during helium burning and can undergo another phase of 
mass loss. Any estimate of mass loss is hampered by the lack 
of non-linear calculations during these late phases of evolution. 
Moreover, as the star evolves toward low effective temperature 
and eventually becomes a red supergiant, it develops a deep 
convective envelope, the opposite of main sequence stars. This 
adds additional difficulty to any derivation of M since convec- 
tion can have an important contribution to the damping or driv- 
ing of the pulsation in this case. 



5. DISCUSSION AND CONCLUSIONS 

Our analysis has shown that even though metal-free stars 
with I2OM0 <M< 500 Mq are radially unstable on the main 
sequence (like their metal-enriched counterparts), the efficiency 
of the e-mechanism is reduced because of the hotter central T 
and more compact structure. Consequently, the e-folding times 
characterizing the growth of amplitudes are much longer After 
burning some hydrogen, the stars stabilize (Table 2). 

In order to know the effects of such oscillations and accu- 
rately estimate the resulting pulsationally driven wind, non- 
linear calculations would be required. We have not done those 
and have only estimated the mass loss using simple arguments 
based upon previous non-linear calculations. Adopting the ap- 
proach of Appenzeller (1970ab) that the mass loss can be de- 
rived from the energy gain rate of the pulsation when its sur- 
face velocity reaches the local velocity of sound at least gives 
a quantitative result which is a reasonable upper limit (Sect. 4). 
For the I2OM0 star we find negligible mass loss during the 
central hydrogen burning phase. The integrated mass loss dur- 
ing central hydrogen burning of about 2.5 Mq for a 200 Mq star 
increases to about 25 Mq for the 500 Mq star - less than 5 % of 
its initial mass. 

We also note that Papaloizou (1973b) disagrees with the re- 
sults of Appenzeller (1970a) and finds an important limitation 
of the pulsation amplitude, which cannot drive any mass loss. 
If this is the case for the solar metallicity stars, it is even more 
so for metal-free stars. Therefore we conclude that mass loss 
during central hydrogen burning remains unimportant for the 
evolution of all the Z=0 stars investigated (Table 3). 

Finally, it is worth noting that the strong instability found at 
the end of core helium burning for M > 300 Mq could have im- 
portant consequences. What exactly happens depends on the re- 
maining time before collapse (only a few 10 kyr for the 300 and 
500 Mq stars investigated here) and, especially, on the strength 
of the mass loss rates achieved in these pulsating giant stars 
with convective envelopes. This is an important issue regarding 
the gamma-ray burst scenario suggested by Fryer et al. (2000) 
and needs further investigation. 

From our estimates of the mass loss rates we conclude that 
primordial very massive stars of < 500Mq (but over IOOMq) 
reach core collapse with helium cores massive enough to en- 
counter the pair creation instability. Single stars probably 



Table 1 
Properties of ZAMS models as a function of mass and metallicity. 



M(Mq) 


z 


Teff (kK) 


R(Rq) 


rc(10«K) 


P(d) 


Td (yr) 


120 





99.2 


4.6 


1.38 


0.064 


1.4x10"* 




2x10-*' 


72.9 


8.2 


0.74 


0.16 


5.6x10^ 




2x10-5 


67.1 


9.6 


0.63 


0.20 


4.2x10^ 




2x10-"* 


62.6 


11.0 


0.55 


0.25 


3.7x10^ 




0.01 


53.3 


15.2 


0.45 


0.34(0.16) 


4.6x10^(2.3x10') 




0.02 


48.9 


18.4 


0.44 


0.36 (0.26) 


1.9x10^(1.4x10-^) 




0.02 (g- scatt.) 


67.1 


13.2 


0.45 


0.34 


7.8 xlO^ 


200 





104.2 


6.0 


1.43 


0.083 


3.7x10^ 




0.02 


44.9 


32.1 


0.46 


0.76 


1.3x10-3 




0.02 (g- scatt.) 


70.5 


17.5 


0.47 


0.45 


5.3x10^ 


300 





107.0 


7.6 


1.47 


0.10 


2.6 xlO^ 




2x10-"* 


68.1 


18.2 


0.60 


0.39 


1.3 xlO^ 




0.02 


39.8 


54.1 


0.47 


1.79 


2.1x10-3 




0.02 (g- scatt.) 


72.5 


22.0 


0.48 


0.55 


4.3x10^ 


500 





109.2 


10.0 


1.51 


0.13 


2.0x10' 



Note. — The pulsation period and e-folding time are given for the fundamental mode. Only for 120 M© with Z > 0.01 the fun- 
damental mode and first overtone (in parentheses) are both unstable. For Z = 0.02, results based on models calculated with pure 
electron scattering opacity {e~ scatt.) are also given. 



Table 2 
Properties of Z = stars during core hydrogen burning. 

M(Mq) TH(Myr) r,t,bie (Myr) YA%) A^,, 



120 


2.35 


0.25 


69 


10 


200 


2.15 


0.52 


58 


69 


300 


1.94 


0.63 


51 


130 


500 


1.74 


0.70 


46 


216 



120es 1.74 0.97 33 1091 



Note. — The table gives the initial mass M, the total hydrogen-burning lifetime of the star tu, the age at which the star becomes 
stable Tstabie, the central hydrogen mass fraction at this time Yc (initial value is 76 %), and the number of e-folds an initial pulsation 
has grown since the ZAMS by then, A^ax = Jq""''" K(t)dt . 



Table 3 
Mass loss estimates of stars at ZAMS. 



M(Mq) 


z 


£pui.s(1048 erg) 


M(M0yr-i) 


AM(Mo) 


120 





0.5 


7.5x10-' 


0.2 


200 





1 


5 xlO-*- 


2.5 


300 





2.5 


1.3 xlO-5 


8.5 


500 





6 


3.5 xlO-5 


24.5 


120 


0.02 (e- scatt.) 


0.2 


1.7 xlO-5 


17 



Note. — The table gives the initial mass, M, metallicity, Z, pulsation energy, fipuis, mass loss rate at the ZAMS, and an estimate of 
the total mass lost during the time the stars are unstable (rstabie, see Table 2) during central hydrogen burning, i.e., /SM = MTstabie- The 
results for 120 M©, Z = 0.02 and pure electron scattering opacity is also given. 
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do not lose their hydrogen envelope for initial masses of < 
300 M0. The nucleosynthetic yields of these stars can be im- 
portant for the chemical evolution of the early universe, as they 
may eject up to 60 M© of ^*'Ni for a « 250Mq star (Heger & 
Woosley 2000). Still heavier stars make black holes (Bond, 
Arnett & Carr 1984). Extrapolating from our results, we specu- 
late that stars of >1OOOM0 might experience such strong mass 
loss that they end their lives as smaller objects. But clearly 
the stability of these massive "first stars" is a subject worthy of 
additional exploration beyond the first linear stability analysis 
presented here. 
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